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A  quantum-kinetic  theory  has  been  formulated  to  describe  the  spatiotemporal  microscopic  carrier  dy¬ 
namics  in  laser  excited  semiconductor  materials  that  accounts  for  the  effects  of  inhomogeneous  excitation 
and  structural  inhomogeneities  due  to  bulk  fillamentation  damage  and  micro/nano  structuring  by  introducing 
space  dependent  carrier  distributions  and  coherency.  Two  alternative  approaches  to  the  problem  have  been 
presented  -  one  is  a  non-trivial  generalization  of  our  previously  derived  Fokker-Planck  equation  in  energy 
space,  and  the  other  is  density  matrix  based  formalism.  Both  approaches  involve  extensive  computational 
effort. 

I.  NON-LOCAL  FOKKER-PLANCK  EQUATION  FOR  LASER-INDUCED  BREAKDOWN  IN 

MATERIALS 

A.  Introduction 

Ultra-short  laser-induced  breakdown  studies  have  revealed  the  possibility  of  self-channeling  ultrafast 
light  in  condensed  media  with  the  strong  involvement  of  single  or  multiphoton  ionization  and  avalanche 
ionization  capable  of  delivering  high  free-electron  densities.  The  full  description  of  the  laser-matter  interac¬ 
tion  process  leading  to  material  modification,  and  the  complex  spatio-temporal  dynamics  of  the  beam  prop¬ 
agation  in  the  media  from  first  principles  comprises  the  nonlinear  Schrodinger  equation  for  the  laser  electric 
field  envelope  coupled  with  an  equation  describing  the  growth  of  the  free-electron  density  n(r,  t)  [1],  [2]. 
Usually  the  equation  for  n(r,  t)  employed  in  the  literature  is  a  simple  rate  equation  not  accounting  for  the 
energy  absorption  process  from  the  laser  field.  The  slowly  varying  [1]  envelope  function  E  of  the  electric 
field  of  the  laser  E  is  defined  by:  E  (r,  t )  =  E  ('A  t)  e +  c.c.  (for  a  pulsed  beam  centered 

on  the  operating  frequency  <  wave  number  /,;q  ) ;  we  choose  the  electromagnetic  radiation  to  be  polarized 
in  the  x  direction  E  =  (E,  0,  0).  The  free  electrons  motion  is  governed  by  the  oscillation  of  the  inhomoge¬ 
neous  optical  field  beating  at  Ql  so  that  an  electron  current  density  with  a  slowly  varying  envelope  induces 
a  variation  of  the  local  population  of  electrons  in  space  and  time  leading  to  a  space-dependent  carrier  distri¬ 
bution.  Since  the  time  development  of  the  avalanche  ionization  process  depends  also  very  strongly  on  the 
values  of  the  field,  it  follows  that  at  any  spatially  non-uniform  distribution  of  the  electromagnetic  radiation 
density,  initial  electron  density  gradients  Vr  could  occur  that  may  give  rise  to  more  or  less  strong  diffusion 
currents.  Therefore,  instead  of  using  a  simple  rate  equation  for  the  electron  density  n{r,  t)  produced  by  the 
laser  field,  we  augment  phenomenologically  the  kinetic  equation  previously  derived  by  us  for  the  electron 
energy  distribution  function  to  account  for  its  spatial  dependence  when  the  laser  field  is  spatially  inhomo¬ 
geneous.  The  rates  in  the  equation  explicitly  include  the  effect  of  laser  field  energy  absoiption  by  the  free 
electrons. 
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B.  Theoretical  model 


A  kinetic  equation  is  derived  for  the  conduction-electron  distribution  function  following  references  [3, 4, 
5].  We  assume  that  the  frequency  of  the  electromagnetic  field  satisfies  the  conditions  1  <C  Q/tp  <C  c/v  (the 
electromagnetic  field  frequency  exceeds  appreciably  the  frequency  of  electron-phonon  collisions),  U  A>  uip, 
£Tp  1  (where  rp  is  the  momentum  relaxation  time,  up  is  the  electron  plasma  frequency  and  E  and  v  are 
the  average  energy  and  velocity  of  the  electron).  These  conditions  are  sufficient  to  guarantee  the  penetration 
of  the  electromagnetic  field  deeply  into  the  semiconductor  and  also  mean  that  the  wavelength  is  much  larger 
than  the  mean-free-path  or  the  De  Broglie  wavelength  of  an  electron,  the  wavelength  of  the  phonons  which 
interact  with  the  electrons  and  allow  for  the  use  of  dipole  approximation  of  the  radiation  field.  They  imply 
that  the  dependences  of  these  processes  on  the  electromagnetic  field  are  local.  Under  the  same  conditions, 
the  magnetic  field  contribution  to  the  interaction  with  the  laser  field  is  negligible  because  it  is  of  the  same 
order  of  magnitude  as  the  electric  quadrupole  contribution.  We  first  consider  an  electron-phonon  system  in 
a  spatially  uniform  sinusoidal  electric  field  with  a  vector  potential  A(t)  =  A  .sin(Ofi)  and  find  the  solution 
of  the  Schrodinger  equation 


ih 


dt 


2m* 


i  2 


p  +  -A(r,t)  ip(r,t ) 


(1) 


Using  the  standard  time-dependent  perturbation  theory  and  the  calculated  electron  wave  function  renor¬ 
malized  by  the  laser  field,  we  calculate  the  electron  transition  rate.  Unlike  the  isotropic  transition  rate 
obtained  in  the  absence  of  the  laser  field,  the  presence  of  linearly  polarized  electromagnetic  radiation  re¬ 
sults  in  an  anisotropic  electronic  transition  rate  characterized  by  field  dependence  via  a  Bessel  function 
term  -  this  effect  is  known  as  free-electron  absorption  of  photons.  For  moderately  high  intensities  of  the 
laser  radiation,  the  Bessel  function  is  expanded  to  second  order.  In  the  diffusive  limit  Hujq  <C  Ef.,  where 
Ep  =  h2k2 /2m*e  and  hwq  are  the  free  electron  kinetic  energy  in  the  conduction  band  with  effective  mass 
m*  and  the  phonon  energy  respectively,  the  electron  distribution  function  is  expanded  in  Taylor  series. 
Thus  the  linearized  Fokker-Plank  type  equation  is  obtained.  The  classical  Joule  heating  of  electrons  from 
the  field  which  is  conspicuously  missing  in  other  quantum-kinetic  approaches  is  included  by  taking  into 
account  the  local  fluctuation  of  the  electron  kinetic  energies.  Source  terms  such  as  stimulated  interband 
electron  transitions  due  to  a  single  or  multiphoton  absoiption,  impact  ionization  due  to  Coulomb  interaction 
between  electrons  and  holes  and  sink  terms  such  as  non-radiative  recombination  due  to  a  phonon-mediated 
interaction  are  included  in  the  equation: 


dfe ( E,t ) 


+  V(E,t) 


dfe(E,t ) 


D(E,t) 


d2fe(E,t) 

HEtf 


dt 


dE 


A  (E,  t)  fe  (E,  t)  +  S  {E,  t)  .  (2) 
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The  rates  in  equation  (2)  arc  given  by  the  following  expressions  [3],  [4]: 


V  (E,  t )  =  Ft  (E)  +  VF  (E)  +  ^  (UL)  E2  +  AT(E)  (fiL)  E2, 

D  (E,  t)  =  Dt  (E)  +  Df  (E)  +  ^ a  (QL)  E 2Eek  +  VT(E)\ j  (fi£)  E2, 


(3) 

(4) 


where 


vT(Ei)  =  F  EIc 


h  ^  I  ix 

~q\ 


Hu) X 

q\ 


Nts  {Ei~Eur^>)  -  (Nt+1) 6  (£r-i 


,  (5) 


^eI^-x^-e^ 


h  ^  i  9- 

~q\ 


6 \Ei- Eu~ h^x+hnrJ+S  \En - Eu- hutiX -  hnL 


+  (n**+i)  (Ek-Ek_q)  [<*  (i 


\Ek~Ek+q  +  fUjq  A  +  ^i)  (J 


“  -Bfc-nf+  ^9  A  “  mL 


(6) 


,  (V) 


q  A 


O'-n  E  IC?aI291  X  iV^-E*-,)2  +  + 


4fi^  9- 

~q\ 

+(NP^  +  l)(Ek  -  Ek_qY  5{E%  -  Eu<r+  fuv^x  +  mL)  +  5{E%  -  EUff+  -  MlL) 


q  a 


k+q 


q  A 


(8) 


The  right-hand  side  of  (1)  contains  the  coefficients  for  thermal  spontaneous  phonon  emission  and  the 
one  for  the  field-induced  phonon  emission. 


Mec)  =  ^  ■£  \c. 


h  ^  '  <1X 

~q  A 


"?/(E£  -  Ej+?-  fc?A)  -  (JV£  +  l)i(EE  -  EE+f  + 


ph 


,  (9) 


.  vr 
—  “ 7; 

2h 


+ 


<?  A 1 


fc-q  q  A 


9  A 


“  Ek+q  +  +  mL)  +  &{E%  -  EA-r+(7+  Tuv-jx  -  hnL) 


q  a 


fc+9 


9  A 


(10) 


In  these  expressions  E  denotes  the  amplitude  of  the  laser  field,  the  indices  T  and  F  denote  the  thermal 


contribution  due  to  phonon  scattering  and  field  induced  contribution  respectively, 
phonon  coupling  and  Ssg  =  I  e 


C 


q  a 


is  the  electron- 


g.E  (t)  j m*fi2  J .  The  expressions  for  the  field  dependent  rates  A F,  VF 
and  DF  reveal  that  in  the  presence  of  electromagnetic  radiation  the  electron-phonon  scattering  can  be  ac¬ 
companied  by  electronic  transitions  through  photon  absorption  and  emission. 
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C.  Generalized  Fokker-Planck  kinetic  equation 


By  explicitly  substituting  the  spatially  dependent  envelope  of  the  laser  E  (r,f)  field  in  equations  (3)- 
(10)  the  space  and  energy  dependent  rates  V  (r,  E.  t)  and  D  Cr,  E,  t)  are  obtained.  Equation  (2)  is  different 
from  the  standard  Fokker-Planck  equation  since  the  derivation  does  not  make  it  a  conservation  equation  of 
the  electron  number  and  current  density  even  in  the  absence  of  the  source  terms.  We  define  the  current  in 
energy  space: 


JE(E,t ) 


V(E,t)  fe  (E,  t) 


d  [D{E,  t)fe  ( E,t )] 
dE 


and  equation  (2)  is  written  as 

‘)  r  <* <*, t)  r  (E, t)+s(E, t). 

For  the  time  being  we  ignore  the  third  term  on  the  left-hand  side  of  equation  (11).  Bearing  in  mind 
that  the  Joule  heating  V  C r,E,t )  and  diffusion  D  (r.  E,  t)  rates  of  equation  (2)  already  include  space 
dependence,  a  generalized  current  J°  (r,  E,  t)  is  defined  reflecting  the  energy  gain  or  loss  by  electrons.  The 
electron  distribution  function  is  also  generalized  to  include  space  dependence 


Ju  (E,r,t)=V(E,r,t)f(E,f,t)~- 


d[D(E,r,t)f(E,r,t)} 


dE 


In  references  [6]  and  [7],  a  theory  of  electron  drift  and  diffusion  was  proposed  allowing  for  a  non¬ 
equilibrium  energy  distribution,  such  as  occurs  in  high  electric  fields.  The  transport  equation  assumes 
a  form  of  a  Fokker-Planck  equation  in  a  four-dimensional  (4D)  energy-position  continuum  instead  of  the 
3D  position  space.  In  this  equation  the  4D  current  density  is  a  sum  of  a  drift  and  a  diffusion  terms  that 
depend  on  the  (E,  r)  point.  The  use  of  this  equation  assumes  that  not  only  the  energy  but  also  the  position 
coordinates  vary  continuously,  therefore  demanding  that  the  electron  energy  be  exchanged  with  the  lattice 
in  small  steps  (e.g.  phonons  of  energy  hujq  <C  E)  and  that  electron  velocity  be  finite. 

Following  [6],  electron  transport  in  an  inhomogeneous  laser  field  involving  continuous  energy  and  coor¬ 
dinate  changes  is  described  by  an  equation  of  continuity  in  the  energy-position  manifold.  So  we  generalize 
equation  (11)  as 

dfe  (E,  r,  t )  dJ°  (E,  f,  t )  d J13  (E,  f,  t) 

dt  +  dE  +  d HP 


A{E,r,t)  +  A(E,r,t)  fe  (E,  f,t)  +  S  (E,  f,  t)  (12) 


by  adding  currents  in  real  space,  where  the  Greek  indices  run  from  1  to  3. 
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Due  to  our  choice  of  laser  field  polarization  in  the  x  direction,  we  get  for  the  currents  in  space: 

d[eE  (r,  t)Dx(E,  r,  t)fe(E,  r,  t)]  d[Dx(E,  r,  t)fe(E ,  r,  t)] 


JX(E,  r,  t )  =  vd(E,  f,  t)fe(E,  r,  t )  -  - 


dE 


dx 


JV{E,  r,t)  =  - 
Jz(E,r,t)  =  - 


d[eE{r,  t)Dx(E,  r,  t)fe{E,  f,  £)]  _  d[Dx{E,r,t)fe(E,r,t)] 
dE  dy 

d[eE(r,  t)Dx{E,  r,  t)fe{E,  r,  i)]  d[Dx(E,  r,  t)fe(E,  r,  i)] 


dE 


dz 


(13) 

(14) 

(15) 


In  equations  (13)— (15),  E  (r.  t)  is  the  laser  field  envelope.  We  also  generalize  the  expression  for  single 
photon  excitation  across  the  bandgap.  The  other  two  terms  included  in  S(E,  f,  t )  -  the  impact  ionization 
and  the  Auger  recombination  -  are  not  affected  directly  by  the  laser  field  in  our  derivation. 


Sabs  OC  ^  \Fk\2 


2Fk/n 


_  (HQl  -El-  El  -EGy  +  4  I  Fk[ 


2  e2E  lL{x,y,z,t) 


m0£lL 


rrip 

m* 


Eg  {Eg  +  Ap) 

2  {Eg  +  2A0/3)_  ' 


(16) 

(17) 


D.  Relation  between  the  rates  in  energy  and  coordinate  space 


Following  closely  references  [8,  9,  10],  we  establish  the  connection  between  the  real  space  drift  ve¬ 
locity  and  the  field  dependent  velocity  in  energy  space  which  we  have  already  derived.  In  this  way  we 
circumvent  the  necessity  to  introduce  an  expression  for  the  energy  dependent  electron  mobility  in  the 
given  material  and  also  use  a  quantum  mechanically  derived  rate.  We  consider  an  electron  with  mo¬ 
mentum  pa  and  energy  E  (pa)  at  time  t  =  0.  We  define  dp  =  e E  •  Vp  where  E  denotes  the  negative 
laser  electric  field.  Then  the  energy  change  after  t  =  0  before  any  collision  occurs  is:  E  (t)  —  Eq  = 
( dFE)po  t  +  1/2  (&fE)  i1  +  . . .  The  probability  for  the  electron  starting  from  pa  not  to  experience  col¬ 
lisions  between  0  and  t  is  c(i)  =  exp  |  —  f  dt' /[r  (p)]|,  where  account  is  taken  of  the  change  in  r  as 


the  electron  momentum  changes  in  time  as  dp  j  dt  =  eE.  Upon  first  order  expansion  of  r  ypa  +  eEtj ,  we 
obtain  c  (t)  =  exp  (-t/r0)  1  +  1  / 2  (t/r0)2  (dFr)pQ  +  O  (E2) ,  where  r0  =  r  (p0). 

The  exponential  gives  the  non-scattering  probability  if  r  is  constant  while  the  electron  is  accelerated  by 
the  field  and  the  second  factor  is  the  lowest-order  correction  due  to  the  change  in  r.  We  obtain  the  averages 
over  a  collision-free  flight 

(f)  =  r0  1  +  (dFr)po  and  ( t 2)  =  2 r02  1  +  3  (. dFr)po  . 

Then  the  average  energy  gained  over  a  collisionless  flight  (of  order  E2)  is: 

{E  (i)  -  Eo)  =  {( dFE)po  +  dF  (^|)po}  r0  so  (E  (i)  -  E0)  =  eE  ■  {ns  (pb)  +  dF  (vsr)ft]  }  r0, 

J  J  E0dF  (vgT  (p))  dSp/\Vg\ 


VF  {Eo)  =  (<§)F  =  eE  •  { Vg  {po)  +  dF  (uflr)ft  }  ,  with  vd  {E0) 


I  I  E0dSp/\Vg 
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where  dSp/\vg\  are  the  number  of  states  in  a  surface  element  dSp  in  momentum  space.  We  thus  finally 
obtain  the  relation  between  the  field  dependent  drift  velocity  and  the  drift  velocity  in  coordinate  space: 
Vf  (Eq)  =  eE  ■  Vd  (Eq)  and  Vdx(E ,  r,  t)  =  The  position-space  diffusion  coefficient  is  related  to 

the  field  dependent  diffusion  coefficient  in  energy  space  by  [8] 


Dx(E,f,t) 


DF(E,f,t) 

(eE)2 


Joule  heating  is  the  main  energy  gain  mechanism,  so  the  electrons  from  low-energy  states  in  the  conduc¬ 
tion  band  are  transferred  to  high-energy  states  by  absorbing  energy  from  the  laser  field.  This  is  illustrated 
in  figure  1  where  the  peak  of  the  electron  energy  distribution  function  exhibits  increased  amplitude  and  is 
shifted  to  the  right  (higher  electron  energy)  in  the  presence  of  the  Joule  heating  term.  Figure  2  shows  the 
dependence  of  the  electron  average  energy  on  the  Joule  heating  term.  It  should  be  noted  that  the  Joule 
heating  term  is  only  one  contribution  to  the  rate  V  ( E,t )  (eq.  3)  which  has  the  physical  meaning  of  energy 
dependent  drift. 


Figure  1.  Comparison  of  the  energy  spectra  of  at 
f/r^  =  —0.5  including  and  excluding  Joule  heating. 


Figure  2.  Average  electron  energy  at  scaled  time  of 
the  laser  pulse  including  and  excluding  Joule  heating. 


E.  Conclusion 

In  our  theoretical  consideration,  we  have  assumed  that  the  wavelength  of  the  laser  radiation  is  longer 
than  the  other  characteristic  lengths  so  in  the  quantum-mechanical  derivation  we  have  ignored  the  spatial 
dependence  of  the  electromagnetic  field.  In  the  initial  quantum-mechanical  calculations  we  have  consid¬ 
ered  an  electron-phonon  system  in  a  spatially  uniform  sinusoidal  electric  field  of  frequency  Ql  (dipole 
approximation).  After  the  derivation  of  the  kinetic  equation  for  the  electron  energy  distribution  function, 


the  spatial  dependence  of  the  laser  field  is  taken  into  account  by  explicitly  including  in  the  rates  of  the  equa¬ 
tion  the  spatially  dependent  envelope  function  that  comprises  the  amplitude  of  the  laser  radiation  electric 
field.  The  left  hand  side  of  the  initially  obtained  kinetic  equation  is  represented  as  a  current  in  energy  space. 
Consequently  the  energy  dependent  current  is  generalized  to  be  energy  and  space  dependent.  The  kinetic 
equation  is  augmented  with  spatially  dependent  currents  and  a  generalized  Fokker-Planck  type  equation 
in  the  energy-position  manifold  is  obtained.  Without  invoking  expressions  for  the  electron  vector  mean- 
free-path  related  to  the  momentum  relaxation  and  the  electron  mobility  of  the  sample,  from  elementary 
arguments  we  have  obtained  the  connection  between  the  field  dependent  drift  velocity  in  energy  space  and 
the  drift  velocity  in  real  space. 

II.  QUANTUM  KINETIC  THEORY  FOR  SPATIO-TEMPORAL  DYNAMICS  OF  CARRIERS  IN 

LASER-EXCITED  MATERIALS 

A.  Introduction 

When  dealing  with  a  typical  semiconductor-based  optoelectronic  device  irradiated  with  laser  light  we 
consider  the  optical  field,  the  created  electron-hole  plasma  (EHP)  and  the  crystal  lattice  of  the  chosen  ma¬ 
terial.  Light  generation,  propagation  and  amplification  determine  the  behavior  of  the  optical  field.  Carrier 
generation  and  recombination,  electrical  conduction  and  diffusion  determine  the  behavior  of  the  formed 
plasma.  The  photon  energy  of  the  laser  field  is  converted  and  conserved  as  kinetic  and  thermal  energy  of 
the  plasma  and  thermal  energy  of  the  lattice  by  creation  and  annihilation  of  phonons.  All  the  described 
processes  take  place  on  different  time  and  space  scales  but  they  should  be  treated  in  a  self  consistent  man¬ 
ner  with  the  appropriate  coupled  equations.  Laser  beam  filamentation,  dynamic  beam  steering,  catastrophic 
optical  damage,  thermal  lensing  leading  to  formation  of  hot  spots  are  cases  requiring  inclusion  of  spa¬ 
tial  variation  in  the  formalism  describing  the  dynamics  of  optically  generated  carriers  interacting  with  the 
phonons  of  the  lattice.  Inhomogeneous  excitation,  bulk  filamentation  laser  damage,  etc.  lead  to  space 
dependent  carrier  distributions. 


B.  Theoretical  model 

When  a  spatially  homogeneous  system  is  excited  by  a  spatially  inhomogeneous  laser  field,  the  dynam¬ 
ical  variables  become  inhomogeneous  and  off  diagonal  density  matrices  have  to  be  introduced.  A  mixed 
momentum  and  real  space  representation  is  most  similar  to  classical  distribution  function  and  is  best  suited 
for  a  comparison  to  semi-classical  kinetics  described  by  Boltzmann  equation.  A  microscopic  density  matrix 
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theory  is  formulated  accounting  for  arbitrary  spatial  inhomogeneities  in  the  excitation  conditions  leading  to 
space-dependent  Boltzmann-Bloch  transport  equations  for  the  description  of  spatio-temporal  dynamics  of 
electrons  and  holes  of  inhomogeneously  excited  materials  such  as  semiconductors  including  the  coherent 
interactions  of  carriers  and  the  laser  light  field  as  well  as  transport  due  to  spatial  gradients  and  electrostatic 
forces.  Only  the  classical  character  of  the  laser  optical  field  is  considered  while  accounting  for  the  quantum 
mechanical  properties  of  the  semiconductor.  Besides  the  interaction  with  the  light  field  other  important 
interactions  occur  in  the  semiconductor  -  Coulomb  interaction  among  the  carriers  giving  rise  to  screening 
and  to  thermalization  of  the  nonequilibrium  carrier  distribution,  as  well  as  interaction  with  phonons  lead¬ 
ing  to  an  energy  exchange  between  the  carriers  and  the  crystal  lattice.  Based  on  typical  length  and  time 
scales  approximations  are  made  with  the  aim  of  obtaining  numerically  tractable  system  of  equations.  We 
follow  the  approach  in  [11]  but  unlike  them  we  treat  all  the  scattering  terms  explicitly  without  resorting 
to  relaxation  time  approximation  [12].  We  also  include  terms  that  lead  to  transitions  between  valence  and 
conduction  band  -  impact  ionization  and  Auger  recombination  [3]. 

We  consider  a  two-band  model  of  an  udoped  semiconductor  such  as  GaAs.  In  the  laser-matter  inter¬ 
action  process  the  physical  variables  that  are  directly  related  to  observables  of  the  system  such  as  optical 
polarizations  and  distribution  functions  are  all  single-particle  quantities  calculated  by  the  density  matrix.  To 
describe  space-dependent  phenomena  a  Wigner  representation  of  the  single-particle  density  matrix  can  be 
used.  In  Wigner  representation  the  space-dependent  distribution  functions  (intraband  density  matrices)  of 
electrons  and  holes  and  polarization  (interband  density  matrix)  are  defined  as 

ph/ 


fe(k,r)  =  E^'r(c£+i  -+  Cjf  i  _),  fh(k,r)  =  E  e*9'r(<%+i -+ <%  i ->  and 
q  q 


wave 


p{k,f)  =  J2eiq-r(d_Ul^cUy), 

Q 

where  ci  and  di (c^  and  d^)  denote  creation  (annihilation)  operators  for  electrons  and  holes  with 
vector,  respectively  and  the  brackets  denote  the  expectation  value  of  these  operators. 

The  single-particleHamiltonian  describing  the  free  carrier  interacting  with  a  classical  light  field  as  well 
as  the  free  phonons  is  given  by: 


h«  =  Y.  etci  Ce+Y,4dtdi+Y.  bi 


E 

k,  q 


Mk)  ■  Et  (t)c+  d+  _  +  jTcv(k)  ■  E-  ( t)d+  c+  _ 

■*  k-\- q  — fc+  7j  q  *  — k-\-Tjq  k-\-^q 


(1) 


where  fi  (  k )  is  the  component  in  the  direction  of  the  laser  field  polarization  of  the  interband  optical  dipole 


matrix  elememt  between  the  electron  state 


c,  k  )  and  hole  state 


v,—k).  The  field  is  represented  by  two 
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counterpropagating  waves  and  the  positive  frequency  component  is  given  by: 

E+  (f,  t)  =  ^  (E+  (r,  t)  eiK*z~iojt  +  E~  (r,  t )  (2) 

and  is  expanded  in  a  Fourier  series 


E+  (r ,  t)  =  J2  Ef  V)  =  E  (*)  e*'?  (3) 

9  9 

In  the  absence  of  an  external  light  field  the  electron  states  are  eigenstates  of  an  ideal  periodic  lattice.  Devi¬ 
ations  from  this  idealized  periodicity  due  to  lattice  vibrations  lead  to  a  coupling  of  the  different  electronic 
states.  This  interaction  is  described  by  the  carrier-phonon  Hamiltonian. 


H?  =  £  C, 

k,k'  ,q 


ct  bffCr  —  ctbic-  dr ,  -+  b+dt dr  +  b+  +  dr ,  - 
L  k+q  q  k  k  9  k+q  k+Q  9  k  k  9  1  k+q 


(4) 


where  the  electron-  phonon  coupling  constant  for  interaction  with  optical  phonons  is: 


In  |2  (  hwLO\ 

1 

1 

r  e2  1 

1  {  2V  ) 

£r  (oo) 

£r  (0) 

Ao  {q2  +  Q2s)_ 

(5) 

er  (oo)  and  er  (0)  are  the  relative  static  and  optical  dielectric  constant,  respectively,  £o  is  the  absolute 
dielectric  constant  of  the  vacuum,  hi ulo  is  the  optical  phonon  energy  and  V  is  the  normalization  volume. 
The  charged  carriers  interact  via  the  Coulomb  potential  V(j  and  the  Hamiltonian  describing  carrier-carrier 
interaction  processes  conserving  the  number  of  particles  per  band  is  given  by: 


nr  =  £  V, 
k,k\q 


1 


—  C-  C, ,  . 

2  fc  k'  k  +1  k~1 


-Cr  _H — d~tdX  dr, ,  -Ar  -  —  c+d+~d 
a  k-0  2  k  k'  k+Q  k~<l  k  -fc'  - 


-k' +cf"k—q 


(6) 


This  carrier-carrier  Hamiltonian  can  be  separated  into  a  mean  field  (Haretree-Fock)  HC^F  part  and  a  re¬ 
maining  part  depending  on  two-particle  correlations  H"nrr.  The  effective  single  particle  Hamiltonian  is 
Heff  =  II o  +  E[ftF.  The  correlation  part  of  the  carrier-carrier  interaction  Hamiltonian  gives  two  phenom¬ 
ena:  scattering  processes  between  the  carriers  and  screening  of  the  bare  Coulomb  interaction. 

The  part  of  the  perturbation  Hamiltonian  that  yields  impact  ionization  and  its  inverse  process,  Auger 
recombination  is  given  by  [3],  [13]: 


H 


cc(c—v) 


y  I ~Me(q)ct  _ct  _d+  ~  cr  +  M*  (q)  ct d  r,cr,  .ci  _ 

L  y  k+q  k'-q  -fc'  k  e  fc  ~k  k  ~Q  k+q j 


+ 


k,k'  ,q 


+  E  [Mh  (?)  dUqdt'-qC-k'dk  +  Mk  ^  dtC-k'dk'-^k+q 


(7) 


k,k'  ,q 


47TC^ 

Me  (q)  =  Vq9q,  where  V'f  — - -  is  the  Coulomb  potential  and  is  the  interband-transition  form  factor. 

V  CQ4/ 
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C.  Generalized  Boltzmann-Bloch  equations 


By  using  Heisenberg’s  equations  of  motion,  the  equations  of  motion  for  the  single -particle  density  ma¬ 
trices  in  Wigner  representation  can  be  derived.  The  effective  single-particle  Hamiltonian  Heff  gives  a 
closed  set  of  equations  for  the  distribution  functions  of  electrons  and  holes  and  by  interband  polarization. 
Being  Wigner  distributions  these  quantities  are  functions  of  space  and  momentum  but  there  is  a  big  dif¬ 
ference  of  time  scales  between  the  momentum  space  and  real  space  dynamics.  Scattering  and  dephasing 
processes  lead  to  fast  relaxation  of  the  microscopic  variables  towards  their  local  quasi-equilibrium  values 
on  a  femtosecond  time-scale  while  the  spatial  transport  happens  on  a  much  slower  time-scale  (10  ps  to  ns). 
Because  of  the  typical  separation  of  time  scales  between  the  k- space  and  r-space  dynamics,  the  influence 
of  spatial  gradients  on  the  k-space  dynamics  is  often  negligible.  However,  some  of  the  scattering  terms 
in  the  equations  of  motion  for  the  distribution  functions  conserve  the  density  of  carriers  and  therefore  the 
density  is  not  influenced  by  the  fast  relaxation  processes  and  its  spatial  transport  cannot  be  neglected.  In 
the  equation  of  motion  for  the  polarization  no  conserved  quantities  exist  and  thus  the  spatial  transport  of 
polarization  is  usually  not  important.  In  principle  the  complete  set  of  equations  required  is  therefore  the 
Maxwell-Boltzmann-Bloch-Poisson  equations  for  the  nonequilibrium  distribution  functions  fa  (jz,  rj,  in¬ 
terband  polarization  p  (jz.  rj ,  electric  potential  <f>  (r),  and  the  laser  field  E  (f.  t),  with  k  and  r  being  the 
two-dimensional  (2D)  vectors  in  reciprocal  (momentum)  space  and  real  space,  respectively. 

Keeping  the  first  order  spatial  derivatives  of  the  distribution  functions  and  neglecting  any  spatial  trans¬ 
port  of  polarization,  the  equations  of  motion  for  electron  and  hole  distribution  functions  are  given  by  the 
generalized  Boltzmann  equations  for  two  band  model  including  the  coherent  interband  transport  contribu¬ 
tions. 


d  f<*(Z  ^u1  dea(fc,r)  dfa(k,r,t ) 
dt  ’  ’  h  dk  dr 


=  jt“(h,n  +  ^r&r u  (8) 


The  lowest  order  contribution  to  the  polarization  is  included,  where  the  spatial  coordinate  enters  only 
as  a  parameter  and  locally  the  dynamics  coincide  with  those  of  the  inhomogeneous  case  and  there  are 
no  transport  effects.  This  lowest  order  picture  is  sufficient  to  describe  pump-probe  experiments  in  which 
filamentation  is  observed. 

^p(k,r,t)  =  -^[£e{k,r,t,)  +  £h(-k,f,t)\p{k,r,t)- 

-  iU(k,  r)[fe(k,  r,  t)  +  fh{-k,  r,  t)  -  1]  +  r)coi  (9) 


e“  (jz,  r^J  =  £a  (jtj  +  (r)  +  <5e"  (jz,  fj  , 
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~eM  1  k  )  = 


7j—  is  the  single  particle  energy,  Sea(k,  r)  =  -J2  fa(k,  r)  VS ^  ^  [v£,_z  ~  v%'- 
zma  £, 


k'—k  9  fc'— A;  /c'— /c-l 

fc' 


is  the  renormalization  of  the  single  particle  earner  energy  due  to  exchange  interaction,  VS  =  — is 

q  £{q,0) 

the  screened  Coulomb  potential.  The  electrostatic  potential  due  to  the  Hartree  terms  in  the  mean  field 


Hamiltonian  satisfies  the  Poisson  equation: 

cP_ 

dr 2  *  v  '  £qV 


k,  r)  -  fh  ( k,  f 


The  generation  rate  in  the  equation  (1)  is  given  as  follows: 

Ra  (jz,  fj  =  i  (k.  r\  p  *  (jz,  f^j  —  Cl  *  (jz.  rj  p  (jz,  r 
where  ft  (jz,  f'j  is  the  renormalized  Rabi  frequency  defined  by: 

hn  (k,  fj  =  n  (fc)  E  {r,t)  +  ^2p(k\  rj  VS-  . 


(10) 


(ID 


(12) 


The  second  term  in  the  above  expression  is  the  internal  field  responsible  for  Coulomb  enhancement. 


D.  Scattering  processes 


Within  a  semiclassical  picture  when  scattering  processes  are  described  in  terms  of  scattering  rates,  the 
scattering  contributions  to  the  equation  of  motion  have  the  structure  of  the  Boltzmann  collision  terms.  The 
electron-phonon  scattering  rates  are  obtained  from  Fermi’s  golden  rule  and  quadratic  or  higher  order  terms 
and  terms  involving  simultaneous  electron-phonon  and  hole-phonon  interaction  have  been  neglected  in  the 
polarization  and  in  the  carrier-phonon  Hamiltonian.  Incoherent  scattering  processes  appear  for  the  first  time 
in  second  order  contributions  [3],  [11],  [14],  [15]. 

Collisional  contributions  in  equations  (1)  and  (2)  lead  to  relaxation  in  the  carrier  distributions  and  decay 
in  the  interband  polarization: 


d_ 

dt 


r  (*’  ’1  o «  =  §tr  (*> ?) „„  +  Ftf°  (*’  f 

v f)  * = e  H,  £-/(*- *  n  - 


(13) 

(14) 

The  first  term  on  the  RHS  of  equation  (3)  depicts  the  scattering  processes  arising  from  the  correlation  part 
of  the  carrier-carrier  Hamiltonian  and  the  third  term  arises  from  the  carrier-phonon  Hamiltonian 


d_ 

dt 


r(k,f)m/LO 


Q 


e.h 


(k,  0)  -  Wt?,kfe,h^  f)(1  "  f£'h ^  ~  ^ 


(15) 
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where  the  scattering  matrices  are  given  by: 


weMaa)  =  2n  ^  le 


k—q,k 


h 


a=e,h  fc/ 


e,h  .  a  _  a  _  e,h 
k-q  k'+q  k'  k 


/"  (j',  j)  1  -  .f°  (k'  +  q,  j 


we,h(LO)  = 

fc-g,fc  ft 


s~ie,h 


(^  +  1)  <5  -  ef  +  ^)  +  — 


a 


e,h 


N"6[£1-q-£T 


(16) 

(17) 


The  Boltzmann  scattering  matrices  -■  for  electrons  and  holes  are  real  quantities  and  the  scattering 

matices  WS  in  the  equation  of  motion  for  the  polarization  are  complex  and  the  real  part  is  related  to 
We.'h_,  .  according  to: 

k—q,k 


Re  (  WS 

k—q,k 


1 

=  -ITT 

2  k—q,k 


1  - 


/e  (fc  -  9,  j 


+  wb-3rik-*n  + 


+  lw 


2  q—k,—k 


1  -  fh  (q-  k,r 


+  w 


—  k,q—kJ 


fh  q-k,f)  (18) 


The  real  part  of  WS  describes  scattering  processes  leading  to  a  dephasing  of  the  polarization  and  the 
imaginary  part  describes  second-order  contributions  to  the  band-gap  renormalization. 

The  carrier-carrier  scattering  rate  in  the  collisional  contribution  to  the  polarization  equation  describing 
the  effect  of  correlations  is  given  by: 


W- 


-p(aa) 


k—q,k 


EEIW 

a=e,h  fc/ 


e.h  i  a  a  e,/i\  v , 

:J  +  £?,  £?,  -  e-  |x 

k-q  k'+q  k'  k 


x  f°  (k',  rj  1  -  /“  +  9,  rj  x  1  -  /e,/l  (ft  -  q,  j 


+ 


+  rh  (k-q,  f)  r(k'  +  q,  r)  [l  -  /“  (k'  + 


q,r 


(19) 


The  carrier- phonon  scattering  rate  in  the  collisional  contribution  to  the  polarization  equation  is  given  by: 


W- 


■p(LO)  _  7r 


k-q,k  ft 


\Cq\2S  ( £tq  -  £t  +  huJi)  {N^£)h  (*  -  j)  +  +  !)[!-  fe,h  (*  -  9,  j]  }  + 

+  |  lOrl2  -5  (VA  -  ef  -  fu+qj  { (Nq  +  1)  rh  (k  -  q,  rj  +Nq[l-  fe’h  (k  -  q,  rj]  }  (20) 


E.  Impact  ionization  and  Auger  recombination 

Since  we  interested  in  the  processes  of  laser  damage  and  filamentation  in  the  semiconductor  material, 
we  include  the  j-/Q  ( k.  r\  terms  that  lead  to  transitions  between  valence  and  conduction  bands,  i.e. 

ut  V  J  eh 

impact  ionization  term  and  Auger  recombination  term  [3] .  Impact  ionization  and  Auger  recombination  are 
second-order  two-particle  Coulomb  scattering  processes  (proportional  to  Coulomb  scattering).  In  a  case 
when  we  have  a  homogeneous  system  (material)  that  is  either  homogeneously  or  inhomogeneously  excited 
the  Coulomb  matrix  elements  depend  on  the  momentum  transfer  q  only,  i.e.  oc  \V$\  . 
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These  contributions  are  derived  using  second-order  perturbation  theory  such  as  Coulomb  scattering  [3] 
which  is  very  important  to  conduction-electron  dynamics  and  the  change  of  electron  density. 


d_ 

dt 


r  -n^t  ?  MV  (2£i  -  + -ie-,1 + £c) x 


1  ~fe(k,r 


1  2 


f  (k  +  q,  r 


l-fn(k-q,f 


+ 


+  iVe//y  NV  (2e| k-^  -£l  +  £\k-2q\  +  Eg)  x 


xf  k,r 


1  -  fe  [k-q,r 
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l-/fc  fc-2g,f 


(21) 
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fe  (k  -  q,  r)  fh  (k 
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2vr 
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Etf-Ec)x 


fe(k,r 


1  2 


1  -  /e  +  q,  r 


fh{ 


k  —  q,r 


(22) 


withpq-«  2  (m*/m0). 

The  semiclassical  generation  rate  for  earner-light  interaction  is  obtained  by  eliminating  the  polarization 
as  independent  variable  [15].  This  is  done  by  solving  the  equation  for  polarization  within  the  adiabatic  and 
Markov  approximation.  For  Gaussian  pulse  with  a  space  dependent  amplitude 


E  (f,  t)  =  El  (r)  exp 


(' t/TL f 


(23) 


the  time  integration  of  the  polarization  equation  (2)  gives 

2 


Rc 


fr  A  m  a/2  (MkEL{f)\  [  , 

[k,  r  J  =  (27r)  '  I  - - - J  tl  exp  ^-2  [t/rL) 


exp 


(TLAwfc)2 


i -r  k,f  - 


fh  (k,  r 

(24) 


where 


Auk  =  (ee  (k,  rj  +  eh  (k,  r)  )  (25) 

is  the  detuning  of  a  given  transition  with  wavevector  k  from  resonance,  lol  being  the  laser  frequency. 

The  time  integration  is  possible  only  under  the  assumption  that  during  the  laser  pulse  the  polarization 
is  not  influenced  by  any  scattering  processes  leading  to  phase-breaking  during  carrier  generation.  The  only 
density  dependence  of  this  rate  is  due  to  phase-space  filling.  While  adiabatic  elimination  of  polarization 
leads  to  a  simple  closure  of  the  total  set  of  equations,  this  set  of  equations  have  a  severe  deficiency,  especially 
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in  the  presence  of  any  kind  of  spatial  inhomogeneity.  That  is  why  in  the  formulated  approach  we  are  keeping 
the  Boltznann-Bloch  transport  equations  for  the  three  distribution  functions  fe  (k,  fj ,  fh  fj ,  p  (k,  fj , 
though  all  transport  terms  involving  explicit  spatial  variation  of  p  (jt,  fj  are  ignored.  For  ultrafast  spatial 
inhomogeneous  processes  such  spatial  terms  in  the  polarization  equations  should  be  important. 


F.  Numerical  procedures  in  progress 

The  microscopic  dynamics  of  the  distribution  unction  and  the  nonlinear  polarization  are  governed  by 
the  equations  of  motion  (1)  and  (2).  The  /.’-resolved  interband  polarization  equations  have  to  be  solved  self- 
consistently  for  all  space  and  time  grid  points.  A  full  kinetic  treatment  of  the  scattering  processes  will  be 
performed  by  using  e.g.  Monte  Carlo  simulations.  Generalized  Monte  Carlo  methods  taking  into  account 
phase  relations  between  different  type  of  carriers  (polarization  effects),  interaction  of  carriers  with  external 
coherent  inhomogeneous  electromagnetic  field  (generation  effects),  and  the  correlation  and  renormalization 
effects  associated  with  carrier-carrier  interaction  can  be  utilized  [16]. 


G.  Conclusion 

A  microscopic  quantum-kinetic  theory  based  on  density  matrix  formalism  [17]  is  formulated  to  describe 
the  processes  of  short  pulse  laser  interaction  with  materials  such  as  semiconductors  accounting  for  arbitrary 
spatial  inhomogeneities  in  the  excitation  conditions  and  other  spatial  phenomena  such  as  filamentation 
of  tightly  focused  femtosecond  laser  pulses,  structural  modification  and  catastrophic  optical  damage.  A 
system  of  Boltzmann-Bloch  transport  equations  are  established  that  include  both  space  and  momentum 
dependence  of  the  electron  and  hole  distribution  functions  and  the  polarization.  Microscopic  electron- 
phonon  and  electron-electron  scattering  terms  as  well  as  scattering  terms  that  lead  to  transitions  between 
valence  and  conduction  bands,  i.e.  impact  ionization  and  recombination  terms,  are  included  explicitly  in 
the  equations.  The  formulated  theory  describes  the  spatio-temporal  dynamics  of  electrons  and  holes  in 
inhomogeneously  excited  materials  including  the  coherent  interactions  of  carriers  and  the  laser  light  field 
as  well  as  transport  due  to  spatial  gradients  and  electrostatic  forces. 
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